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We derive the optimal N ^ M phase-covariant quantum cloning for equatorial states in dimension 
d with M — kd+N, k integer. The cloning maps are optimal for both global and single-qudit fidelity. 
The map is achieved by an "economical" cloning machine, which works without ancilla. 
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I. INTRODUCTION 

In quantum information the study of optimal cloning machines is a focus of interest since, by definition, cloningis 
synonymous of multiplexing quantum information, which has limitations in principle by the no-cloning theorem 0, 12| ■ 
In the large variety of proposals for optimal doners, the fidelity of the machine depends on the choice of input states, 
with the machine often working in a covariant way, producing "rotated" clones from rotated inputs. In particular, the 
case of SU(d) covariance corresponds to universal cloning with equal fidelity for all unitarily connected states, 

e. g. for all pure states. Clearly, by taking smaller input set of states the cloning performance can be improved, e. g. 
for smaller covariance groups. Also in connection with the eavesdropping strategies in BB84 quantum cryptography 
0, thephase-covariant cloning of equatorial states has been extensively studied for qubits ^^nd more generally 

qudits [3, the latter also with the motivation of understanding which features are peculiar of dimension two. 

In this paper we will consider multi-phase covariant cloning transformations on "equatorial" states 

mq^,})) = ^(1 0) + I 1) + I 2) + ... + e^f^-^ \d-l)), (1) 

where the (f>j^s are independent phases in the interval [0, 27r). An issue which recently has attracted interest in the 
literature is the possibility of achieving the cloning without the need of an ancilla — a so-called "economical" cloning 
As we will see in the following, the multi-phase covariant cloning machines are indeed economical for M = kd+N 
output copies, k integer. 

The paper is organized as follows. In Section^ after introducing the notations and the basic definitions taken 
from Refs. '^^ describe the general approach to covariant cloning maps of Ref. ^^^^ apply it to the case 

of M phase-covariant doner. In Section UTTl we give a brief formalization of the concept of "economical maps" 

by means of the Stinespring representation theorem for completely positive maps. In Section IIVI we explicitly find 
the optimal 1 M doner for both single-qudit fidelity and global fidelities. In Section we generalize all previous 
results to the case A^ M. Section IVll concludes the paper with a comparison of fidelities in the various cases. 



II. PHASE-COVARIANT CLONING 



We want to derive the optimal N M cloning transformations C that are covariant under the group of rotations 
of all the d — 1 independent phases {(pj}, (pj G [0, 2tt), 

t^({'^,}) = |0)(0|+xilj)0>*, (2) 
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where {| 0), | 1), | 2)...| d — 1)} represents a basis for the d-dimensional Hilbert space Ti, of the system of a single copy. 
We win restrict the study of such maps to the set of the TV-fold tensor product of generalized equatorial pure states 

u{{<i^jm^o) = m<Pj})), (3) 

with given in Eq. Here | t/jo) is the equatorial superposition 

i 

The choice 0o = is not restrictive, since an overall phase is negligible. 

As argued in Ref. we consider cloning maps for which the A^-copy input state and the A/-copy output state 
are both supported on the symmetric subspaces Ti.®^ and 7i®*^, respectively. We choose orthonormal basis in the 
symmetric subspace of the form 

= |no,ni,7i2,...nd_i) = -^y Pj^^l 00 . .^11. . . . rf - 1 . . . - 1 ), (5) 

where P^^'^ denotes the permutation operator of N qubits, ng is the number of qudits in state | 0), ni in state | 1), 

and so on, with the constraint X^iLo^ ~ ^ ^'^^ ^^'^ input state, and, analogously, for the output state. In the whole 
paper we will consistently use letters n's for input and rn's for output. The covariance condition for the cloning 
transformation C under the group of multi-phases rotations reads 

As proved in Ref. the covariance condition can be conveniently studied in terms of the positive operator on 

i?=(C®J)(|l))((l|), (7) 
where 2 is the identity map and is the non normalized maximally entangled vector in 7i^^ ® TL®^ 

|1)) = EIM)®IM)- (8) 

The correspondence C R between completely positive maps and positive operators is one-to-one, and can be 
inverted as follows 

C{0) = Tr^«iv [t^^M ® O^) R , (9) 

where denotes the transposition of the operator O with respect to the orthonormal basis in Eq. ^ . Notice that 
for the state jV'o) of Eq. Q one has (|'0o)(''/'ol®^)^ = l'0o)(V'ol'^^ since |'!/'o)j by construction, has all real coefficients 
with respect to the basis in Eq. (jHJ. The trace- preservation condition for C reads 



Tr.^®Af [i?] = (10) 

Following Ref. [T^ , the covariance property ^ rewrites as a commutation relation 

[i?,[/({,/>,})«^®C/*({^,})n =0, (11) 

where the complex conjugated U* of U is defined as the operator having as matrix elements the complex-conjugated 
matrix elements of U with respect to the same orthonormal basis in Eq. (jH)). Eq. (|11|) in turn implies by Schur 
Lemma a block-form for R 

R=^R{m,}, (12) 
{ruj } 

where each set of values {raj} identifies a unique class of equivalent irreducible representations of t^({0j})®*^ ® 
U*{{(t>j})'^^ ■ The equivalent representations within each class can be conveniently written as 

< I Too + no,m-i + ni,nT-2 + ?^2, ■•■TOd-i + no, ni, 712, ...nd-i) f , (13) 

1. J {m] 
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with X]i=o ^ '^j=o "^i = M N. The multi-index {n^} runs over all orthonormal vectors of the basis for 

nf^ used in Eq. ©. With this notation, Eq. ^ becomes 

^=E E ^\7un:^}\{^y} + W^){{^,} + {<}\^\{^){{<^^^ (14) 

{mj} {«;},{""} 

In the following, in order to evaluate the optimality of the map, we will use as figures of merit the single-qudit 
fidelity 

Tr [(|Vo)(V'o| ® 1'-'-' ® |V'o)(^o|^'^) R] (15) 

and the global fidelity 

Tr[(|V'o)(V'o|^*^+^)i?]. (16) 

Notice that in deriving the last two equations we used the covariance property ^ of the input states, the reconstruction 
formula iPJ, the commutation property Hll|l . and the cyclic invariance of the trace. Since each single contribution to 
the single-qudit fidelity H15() and to the global fidelity 1)16(1 is positive versus the indeces {n^} and {n"}, as we will 

show in the following (see Eqs. (|S7jl . and ((1S|| . 0), the block R[m.y must have positive elements — 0^ 

with the off-diagonal ones as large as possible, i.e. '''\™'}{n"} ~ \l ^{n^ }^{n' } \l ^{n;'} {«"} ' E3l- This is equivalent to 
say that the blocks constituting the operator R are actually rank-one blocks, namely 

^{m,} oc |r{„^})(r{„^.} I, (17) 

with 

= E ^{M I {"^.-i + {"^}) ® (18) 
{".} 

and, separately imposing condition (|10|l over every block, Tr.^»Af [| ''{mj})(''{7nj} I] = we get the final form 

for i? 

R= ^ P{m,}\r{m,}){r{m,}l V {m,}) ^ ^ \ [mj] + {n,}) ® \{ni}) , (19) 

{ruj} {m} 

where P{mj} S'l'e free parameters satisfying P{mj} ^ X^Pfmj} ~ 1 order to preserve normalization and 
positivity of R. This means that i? is a convex combination of orthogonal rank-one blocks. 

In Sections II VI and Ivl we will explicitly optimize the map starting from the R operator in Eq. H19() . 



III. ECONOMICAL MAPS 



Let At be a completely positive, trace-preserving map from states on Ti. to states on /C. The Stinespring represen- 
tation Theorem [la | says that for every completely positive trace-preserving map it is possible to find an auxiliary 
quantum system with Hilbert space C and an isometry V from Ti to JC® C, V'^V = 1^^, such that 

M{p)^i:Yc[VpV^]. (20) 

Starting from Eq. (|20|l , it is always possible to construct a unitary interaction U realizing M [T^ ITsI : 

M{p)^Tvc[U{p®\a){a\)U^], (21) 

where \a) is a fixed pure state of a second auxiliary quantum system, say £', such that Ti® C ~ K,® C. The Hilbert 
spaces £ and £' are generally different, and actually play different physical roles. 

We define a trace-preserving completely positive map Jvi to be economical if and only if it admits a unitary form 
U as 

M{p)^U{p®\a){a\)U\ (22) 
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namely, if and only if the map can be physically realized without discarding resources. We can simply prove that the 
only maps admitting an economical unitary implementation U as in Eq. (|22|l are precisely those for which 



M{p) = VpV^ 



(23) 



for an isometry V, V^V = 1. In fact, U{lT-c(^\a)) is an isometry from H to 1C®C, since (l^^® (a|)C/^J7(l-H® |a)) = 
On the other hand, from Eq. H23|l via Gram-Schmidt one can extend any isometry V from 7i to /C (8) £ to a unitary 
U on the same output space, and write it in the form V = C/(1-h ® \o)) for unit vector \a) € with H.® C = JC® C 
For a detailed discussion about the explicit construction procedures, sec Ref. fl^ . 

Allowing classical resources "for free" , a map should be defined as economical also in the case in which it admits a 
random- economical realization as 



Al(p)==^p,[/,(p^|a)(a|)C/, 



(24) 



where pi > 0, J^iPi — 1- Using the same fixed ancilla state \a) for all indeces i is not a loss of generality, since in 
constructing the operators J/i's there is always freedom in the choice of the vector \a). According to this more general 
definition, all economical maps can always be written as a randomization of Eq. (|23|l : 



(25) 



IV. OPTIMAL 1 ^ M CLONING 



The fidelity of the reduced output state Tym-i [C{\ip{{(f>j})) {^p{{(f)j})\)] with respect to the input state 
is given by 



Tr 



h/>o)(V'o|®l'^"'® |V'o)(V'o|) R ■ 
In the case of 1 ^ M cloning, the R operator in Eq. (|19|) has the following structure: 



R= 0P{m,}k{m,})(r{™,}|, 



(26) 



(27) 



with rrij ~ M — 1, whence the form of the summands: 



Tr 



I '0o)(V'o I ® 1*^ ^ ® I V'o)(V'o \ j\mo + l,TOi,m2, ...m^-i) (mo, mi + l,m2, ...m^-i | ® | 0)(1 | 
1 (M-1)! /(mo + l)!mi!....md_i! mo!(mi + l)!....md-i! 1 



■y/{mo + l)(mi + 1). 



(28) 



d2 mo!mi!...md-i! V Ml Ml Md? 

The final contribution to the partial fidelity due to the set of equivalent representations labeled by {m^} is 



^{m,} - Tr 



V'0)(V'0|®]1*'"'® |V'0)(V'0|) R{.n,} 



d2 



1 

M 



^ + !)('"& + 1) 



(29) 



The projector | ''{mj}) (''{mj} I that contributes most to the fidelity is the one that maximizes the quantity 
{mi + l){mk + 1), with the constraint X]f=o^ mi ~ M — I. In the case 



A'/ = dfc + l, /ceN, 



(30) 



the optimization gives the simple result mi = k for all i. The 1 — > {kd + 1) optimal phase-covariant cloning machine 
is then completely described by the rank-one positive operator R = \i'{k}){'''{k}\- The Kraus form of the optimal map 
is then reconstructed as 



C{p) = Tr^ (l^»M ®p^^R = VpV\ 



(31) 
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where V -.Ti. 7i®^^ is the isometry, i.e. V^V = 1, defined as follows: 

V\i) = |mo = fc, mi = k, . . . ,mi = fc + 1, . . .). (32) 

The fact that the Kraus operator describing the map is isometrical — a consequence of R being rank-one — 
automatically guarantees that no additional ancillae (other than the M — 1 blank states) are needed in order to 
unitarily realize the cloning transformation C. 

From Eq. 1)29(1 one obtains the single-qudit fidelity of our multi-phase-covariant economical cloning machine from 
one to M = kd+ 1 copies 

Notice that the above result, in the limit Af — > cx), is consistent with the fidelity of optimal phase estimation on a 
single qudit as worked out in Ref. |lO|| . 

An important remark that remains to be stressed is that the value of {mj} maximizing the single-qudit fidelity 129() 
maximizes the total fidelity as well. In fact, the total fidelity of the {mj}-th block is given by 

[lV'o)(V'o|®('''+'' I . . . ,m, -f- 1, . . .)(. . . ,m, + 1, . . . I ® |z)(j|] = E 1(^0®^*''"'^ I 

ij 

- mo! . . . (m. + 1)! . . . " d^+i ^ Wo; • • • ; + 1; • • -y 

(34) 

where, in the last line, we used the standard notation for multinomial coefficients with the implicit constraints 
= ~ 1- III order to maximize the value of the multinomial coefficient, the vector {nij} has to be as "flat" as 
possible, namely, with all entries as close as possible to each other. The situation in which the solution is unique and 
given by m^ = fc for all i is the same as in Eq. (|3Q|I . This means that the single-qudit fidelity optimization procedure 
provides the same result as the total fidelity optimization, and the map written in Eqs. ((31|l and ((32|l is optimal in 
both approaches. 

V. OPTIMAL N -> M CLONING 

In the general case of arbitrary values for N and M, the single-qudit fidelity is obtained by summing up contributions 
of the form 



Tr 



|^o)(^o| ® ® |^o)(V'o|®'^ \{mj} + {n'}){{m,} + {<}| (g> |K})({<}|1 (35) 



because of the block-form H19|l of the R operator. Before getting into the explicit calculation for the partial fidelity, it 
is possible to somehow simplify the problem by noticing that the presence of \ipo){ipo \ ® l®'*^"^) in Eq. 1(35(1 restricts 
the evaluation of the fidelity to those blocks for which the M-particles states differ at most for a single-particle state. 

The diagonal contributions to F^mj} (as before, the single-qudit fidelity calculated only for the {mj}-th block) are 
then proportional to (apart from the probability p{m^ } ) 



Tr 



\^o){^o\ ^ ^ Ho){'ipo\®^ \ ... ,m, + n, + 1, ...){... + n, + 1, .. . \ ^ \ ... ,n, + 1, ...){... ,n, + 1, . 

1 iV! 1 (?no+?io)!...(m, + ni + l)!... All 

~ d'^ no\...{n^ + iy.... d Ml {niQ + no)l ... {rrii + + 1)1 .. . 

1 A^! 



d^+i no!...(ni + l)!...' 

(36) 

where for sake of simplicity in the last equation the notation was slightly modified with Ui = N — I while, again. 
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M — N . The ofT-diagonal terms are 



Tr 



® ® |Vo)(i^o|®'^ I . . . ,m, + n, + 1, . . .)(. . . , mj + nj + 1, . . . I ® I . . . ,n. + 1, . . .)(. . . , rij + 1, . . . I 



Tr 



li^oXV'ol ® U®(*^-i) I . . . , m, + + 1, . . .)(. . . , mj + + 1, . . . I 



no!...(n,; + l)!... . . . {n, + l)\ 



no! 



(rii + l)!... no!...(nj + 1)!.. 



1 AT! 




/ (mo + no)! . . . {irn + n^ + 1)! . . . 


d^+i no!...7i,!...nj!...\/ 




M! 



(nio + no)! . . . (m^ + n^ + 1)! . 



(M-1)! 



A/! 



(mo + no)! . . . (m^ + n^)! . . . (m^ + n^)! , 



TV! 



Afd^+i no!...nj!...n,!. 



' {rrii + n; + l)(mj + n^ + 1) 

KTT)K"-rT) ■ 



At the end, the single-qudit fideUty is the sum of contributions of the form 



E 



A^! 



?7,o! . . . {n, + 1)! . 



-y 



Nl 



UqI . . .nj , 



{nil + rii + l)(mj + rij + 1) 
(n, + l)(nj- + 1) 



(37) 



(38) 



As done before for the 1 — > M cloning, in order to find the block of R realizing the optimal map, we have to maximize 
the off-diagonal quantity 



(m.j + 111 + l)("ij + nj + 1) 
(^^~TT)(^7~TT) 



(39) 



with the constraints n^ = iV — 1 and mj = Af — N . The maximization of fidelity corresponds to maximize the 
quantity in Eq. (|39() versus the variables m^'s. Since the variables n^'s are summed up in Eq. (|38() . then the fidelity is 
invariant under their permutation. Therefore the evaluation of the maximum of the quantity H39(l resorts to maximize 
it for equal nj's, whence also all m^'s will be equal, m^ = m^ = [M — N)/d, Vz. Generally, in this way one obtains a 
non integer value of m* , while the maximum for integer m^ 's is very degenerate, since the maximum will be obtained 
for unequal m^'s in place of a common fractional value. This leads to many blocks for R contributing in the same 
way to the optimal map, which makes the evaluation very complicate. On the other hand, the evaluation simplifies 
greatly when the maximum is achieved for integer m^ = fc, and this corresponds to the following relation between N 
and M 

M = kd + N. (40) 
Hence, the optimal phase-covariant N —>■ (N + kd) cloning map is described by the rank-one operator 

(41) 



where 



!'"{'=}) E l^"' 



no, 



R = \r{k]){r{k}\, 
-f 111, . . .)M\nQ, .. .,ni. . .)n, E^^J ^' 



and its singie-qudit fidelity is given by 
Fc{N,M = kd + N) = U-^J2i: 



Nl 



d A/d^+i ^ ^ no\...ni^ 



' {n, + k + l){nj + k + 1) 
(n, + l)(?ij + 1) 



E^ 



(42) 



N-1. (43) 



Notice that being R rank-one, the optimal map derived here is again described by only one isometric Kraus operator 



C{p' 



Tr.^«jv 



(U„«„ ® (p®^)^) \r[,}){r{,}\\=Vp^^V\ 



(44) 
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where the isometry V acts as follows: 

V\no, rii, . . . , rii, . . .) n = {hq + k,ni + . . . ,ni + k, . . .)m- (45) 

Similarly to the case 1 (kd + 1), the fact that the optimal N {N + kd) cloning map is isometrical implies that 
no additional ancilla is needed to unitarily realize the map other than the M ^ N blank copies, and (|45|l is again an 
economical cloning machine. 

As in Section llVI it is possible to prove that the value of {mj} maximizing the single-qudit fidelity maximizes the 
total fidelity as well 

J2 E Tr [mi^of^''^''^ + {n'^}){{m,} + {n'^}\ ® \{n'^}){{n'^}\ 

{nij} {n'.},{n'.'} 

= X! X! k'/'^^^^^^'l'^o + no, mi + Til, . . .)|7io,"i, • • 



{nij} {rij} 

1 MT^ / M \ N 

1^ P{m,} 2^ 



(46) 



dM+N + m; • • 7 \no; ni; ■ ■ 

{rnj} {Uj} 

with the usual constraints '^i ^ ^ ^^'^ X]j = M — N implicit in the multinomial notation. Following the 
argument of the previous Section, it is clear that the map in Eqs. H44|l and (|45|l maximizing the single-qudit fidelity 
(|38|l . also maximizes the global fidehty 

As already noticed in the previous Section for N — 1, the fidelity (|43|l . in the limit of an infinite number of output 
copies, namely k — + oo, takes the form (in the limit, M sa kd) 

The above expression coincides with the fidelity of optimal multi-phase estimation on equatorial qudits derived in 
Ref. [13. 



VI. CONCLUSIONS 



We have addressed the problem of optimal phase-covariant cloning with multiple phases for qudits, with arbitrary 
number of input copies N and output copies M. The optimization greatly simplifies for values of M and N related 
a.s M = kd + N , with k integer. The cloning maps are optimal for both global and single-qudit fidelity. The map is 
achieved by an economical cloning machine, which works without ancilla. We have evaluated the asymptotic behaviour 
of the fidelity for large M, and recovered the fidelity of optimal multi-phase estimation In Figs. [3 and 13 

it is possible to compare the single-qudit fidelities of multi-phase covariant and universal covariant I M cloning 
machines. Increasing M the fidelities tend to the limit of optimal phase estimation and state estimation fidelity, 
respectively. Increasing the dimension d of the Hilbert space, the quality of the clones gets worse, as plotted in Fig. 21 
Actually, for fixed N and M, the single-qudit fidelity of the doner goes to zero with d~^, as turns out from Eq. H43|l . 
On the other hand, for fixed M, the fidelity saturates to one as N gets close to A/, in both multi-phase and universal 
covariant frameworks, as one can see in Figs. |31E1 and0 As general remark, notice that, increasing the dimension of 
the input system for fixed N and M, the fidelities of multi-phase and universal doners become closer to each other. 
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